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ABSTRACT 

In this work we present a study of the strong optical collisional emission lines of Ne 
and Ar in an heterogeneous sample of ionized gaseous nebulae for which it is possi- 
ble to derive directly the electron temperature and hence the chemical abundances of 
neon and argon. We calculate using a grid of photoionization models new ionization 
correction factors for these two elements and we study the behaviour of Ne/O and 
Ar/O abundance ratios with metallicity. We find a constant value for Ne/O, while 
there seems to be some evidence for the existence of negative radial gradients of Ar/O 
over the disks of some nearby spirals. We study the relation between the intensities 
of the emission lines of [Nem] at 3869 A and [Om] at 4959 A and 5007 A. This 
relation can be used in empirical calibrations and diagnostic ratios extending their 
applicability to bluer wavelengths and therefore to samples of objects at higher red- 
shifts. Finally, we propose a new diagnostic using [Oil], [Nem] and H5 emission lines 
to derive metallicities for galaxies at high z. 

Key words: ISM: abundances - HII regions: abundances 



1 INTRODUCTION 

The chemical history of the Universe can be investigated 
by studying the behaviour of abundance ratios of different 
chemical species as a function of metallicity, which is the 
main indicator of the chemical evolution of a galaxy. If two 
elements are produced by stars of the same mass range, they 
will appear simultaneously in the interstellar medium (ISM) 
and hence their relative abundance will be constant. But if 
they are produced by stars of different mass ranges, they will 
be ejected into the ISM in different time scales. The chem- 
ical abundances of the elements heavier than hydrogen can 
be studied by measuring the fluxes of the absorption and 
emission lines in the spectra of stars and galaxies. Unluck- 
ily, only the collisional emission lines emitted by the ionized 
gas surrounding massive star clusters are detectable in most 
of galaxies. Since the brightest emission lines in the optical 
part of the spectrum are emitted by oxygen, this element 



* Post-Doc fellow of Ministerio de Education y Ciencia, Spain; 
enrique.perez@ast.obs-mip.fr 

t PhD fellow of Ministerio de Education y Ciencia, Spain 
t On sabbatical leave at IoA, Cambridge 



has been taken as the main tracer of metallicity for these 
objects. Nevertheless, the depletion of some of the most im- 
portant elements, including oxygen, onto dust grains, whose 
composition is difficult to ascertain, makes the determina- 
tion of these abundances more uncertain. This is not the 
case for the elements that occupy the last group in the pe- 
riodic table, that have an electronic configuration with the 
outer shell completely filled and they seldom associate to 
other elements and do not constitute part of the dust grains 
in the ISM. Therefore, although the presence of these grains 
can affect the determination of metallicity in other ways, 
affecting the photoionization equilibrium of the gas, the un- 
certainty due to depletion factors has not to be considered in 
the determination of the chemical abundances of the noble 
gases. 

Neon and argon are products of the late stages in the 
evolution of massive stars. Neon is produced by carbon burn- 
ing and is expected to track oxygen abundances very closely. 
The measurement of neon abundances in extragalactic HII 
regions (Garnett, 2002) and planetary nebulae (Henry, 1989) 
confirm this trend, despite the uncertainties in the deriva- 
tion of the ionization correction factor of Ne. On the other 
hand, argon, like sulphur, is produced by oxygen burning 
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and, again, it is expected to track O abundances. Neverthe- 
less, as it is also the case for sulphur, there is some evidence 
of decreasing values of Ar/O for higher metallicities (Gar- 
nett, 2002). The same trend has been observed for sulphur in 
halo metal poor stars (Israelian & Rebolo, 2001) and extra- 
galactic HII regions (Diaz et al., 1991; Garnett, 2002; Perez- 
Montero et al., 2006). This problem, perhaps, is related to 
the proximity of the production site of these elements to the 
stellar core and the yields would be sensitive to the condi- 
tions during the supernova explosion (Weaver & Woosley, 
1993). 

The emission lines of Ne and Ar are not as intense as 
some of the other strong lines in the optical spectrum but 
there is a growing number of HII regions for which there are 
measurements with good signal-to-noise ratio. This is the 
case of [Nem] at 3869 A, whose blue wavelength, makes it 
observable in the optical spectrum of bright objects, even 
at high redshifts, and the emission line of [Arm] at 7136 A. 
The emission lines of these elements in the IR will supply 
a great deal of worthy additional information but, at the 
moment, the objects with observations of the fine-structure 
lines of [Ne n], [Ne ill], [Ar n], [Ar ill] are still very scarce. 
The derivation of the chemical abundances of these elements 
in HII regions can be calculated by means of the previous 
determination of some electron collisional temperature and 
the previous knowledge of the ionization structure of the 
clement whose abundance is required. The results of pho- 
toionization models point to a similar ionization structure 
of Ne and O on the one hand and S and Ar on the other (e.g. 
Perez-Montero & Diaz, 2007). These facts make these lines 
suitable to be used as substitutes of the bright emission lines 
of oxygen and sulphur in empirical indicators of metallicity 
(e.g. Nagao et al., 2006) or diagnostic ratios to distinguish 
starbursts galaxies from active galactic nuclei (e.g. Rola et 
al., 1997). 

In this paper, we investigate the element abundances 
of neon, argon and oxygen for a sample of objects that is 
described in Section 2. In Section 3, we derive the physical 
conditions of these objects, including the calculation of the 
electron temperatures involved in the derivation of Ne and 
Ar abundances. In Section 4, we describe the grid of pho- 
toionization models that we have used to derive a new set of 
ionization correction factors for these two elements. These 
new ICFs are described in the next Section, along with the 
discussion of the behaviour of Nc/O and Ar/O ratios and 
the use of the brightest emission lines of these two elements 
as empirical calibrators of metallicity and diagnostic ratios. 
Finally, we summarize our results and we present our con- 
clusions. 



2 DESCRIPTION OF THE SAMPLE 

The sample includes emission-line objects observed in the 
optical part of the spectrum with the detection of, at least, 
one auroral line with high signal-to-noise ratio in order to 
derive electron temperatures directly and, hence, ionic chem- 
ical abundances with less uncertainty. We have compiled for 
this sample the strong lines of [On] at 3727 A and [Oiii] at 
4959 A and 5007 A in order to derive oxygen abundances. 
We have compiled as well the [Nem] emission line at 3869 
A, the [Arm] emission line at 7136 A and the [Anv] at 4740 



Table 1. Bibliographic references for the emission line fluxes of 
the compiled sample 



Reference 


Ubject type 


• 

Objects 


Bresolin et al., 2004 


M51 GEHRs 


10 


Bresolin et al., 2005 


GEHRs 


32 


Bresolin, 2007 


M101 GEHRs 


3 


Crockett et al., 2006 


M33 GEHRs 


6 


Garnett et al., 2004 


M51 GEHRs 


2 


Guseva et al., 2003a 


SBS 1129+576 


2 


Guseva et al., 2003b 


HS 1442+650 


2 


Guseva et al., 2003c 


SBS 1415+437 


2 


Hagele et al, 2006 


HUG 


3 


Hagele et al., 2007 


HUG 


7 


Izotov et al., 1997 


IZwl8 


2 


Izotov & Thuan, 1998 


SBS 0335-018 


1 


Izotov ct al., 1999 


HUG 


3 


Izotov ct al., 2001 


HUG 


2 


Izotov ct al., 2004 


HUG 


3 


Izotov & Thuan, 2004 


HUG 


33 


Izotov ct al. 


SDSS galaxies 


309 


Kniazev et al. 


SDSS galaxies 


11 


Lee et al., 2004 


KISS galaxies 


13 


Melbourne et al., 2004 


KISS galaxies 


12 


Perez-Montero & Diaz, 2005 


All 


361 


Van Zee, 2000 


UGCA92 


2 


Vermeij et al., 2002 


DRH 


5 



a GEHR denotes Giant Extragalactic HII Regions, HUG, HII 
Galaxies and DRH, Diffuse HII Regions 



A to measure the respective ionic abundances of neon and 
argon. The compilation includes the objects used in Perez- 
Montero & Dfaz (2005), including HII regions in our Galaxy 
and the Magellanic Clouds, Giant Extragalactic HII regions 
and HII Galaxies, with the addition of new objects which 
are listed in Table 1 

This list includes 12 HII Galaxies from the Sloan Digi- 
tal Sky Survey (SDSS§) with very low metallicities identified 
by Kniazev et al. (2003) and whose spectra have been re- 
analysed to obtain the required line emission data. Among 
these galaxies, 11 do not have observations of the [On] emis- 
sion line at 3727 A, which is the case as well for 183 SDSS 
galaxies of the sample compiled by Izotov et al. (2006). The 
compilation has therefore 633 HII Galaxies, 176 GEHRs and 
44 HII regions of the Galaxy and the Magellanic Clouds with 
a good determination of oxygen, neon and argon ionic abun- 
dances for those ions that have been observed. 

The emission line intensities have been taken directly 
from the literature once reddening corrected. The aperture 
effects of all the compiled observations are negligible due to 
the compact nature of the observed objects and the high ex- 
citation state of the more important involved emission lines. 



§ The SDSS site is at www.sdss.org 
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Table 2. Bibliographic references for the collision strengths of the 
forbidden lines for each ion. 



Ion 


Reference 


On 


Pradhan, 1976 


Om 


Lennon & Burke, 1994 


Nem 


Butler & Zeippen, 1994 


Sn 


Ramsbottom, Bell & Stafford, 1996 


Sm 


Tayal & Gupta, 1999 


Ann 


Galavis, Mendoza & Zeippen, 1995 


Ariv 


Zeippen, Le Bourlot & Butler, 1987 



3 PHYSICAL CONDITIONS 

3.1 Electron density and temperatures 

In order to derive oxygen, neon and argon abundances for 
the ions presenting the corresponding strong emission line 
we have selected objects allowing a direct determination of 
the electron temperature. All physical conditions have been 
derived from the appropriate ratios of emission lines and 
using fittings to the values obtained with the task TEM- 
DEN, in the case of electron densities and temperatures, 
and IONIC, in the case of ionic abundances of the software 
package IRAF ^ and which are based on the five-level sta- 
tistical equilibrium model (De Robertis, Dufour & Hunt, 
1987; Shaw & Dufour, 1995). The references for the colli- 
sion strengths of the ions involved in our calculations are 
listed in Table 2. 

Next, following the same procedure described in Perez- 
Montero &: Dfaz (2003), we have used for the calculation of 
the chemical abundance the electron temperature associated 
with the zone where each ionic species stays. To do this we 
have applied different expressions for the relations between 
the electron temperature of each zone. 

Electron densities have been calculated for a subsam- 
ple of 361 objects from the ratio of [Sn] emission lines 
I(6717A)/I(673lA). It is representative of the low excita- 
tion zone of the ionized gas, and it is therefore used for the 
calculation of electron temperatures in this zone which de- 
pend on density. For the remaining the objects, for which no 
data about [Sn] lines exist or for which they give non valid 
values, we have considered a mean density of 50 particles 
per cm 3 , according to the low density values found in the 
rest of HII regions. 

The electron temperature of 2+ , representative 
of the high excitation zone of the ionized gas, has 
been derived from the ratio of [Om] emission lines 
(I(4959A)+I(5007))/I(4363A) for a subsample of 771 ob- 
jects, while for the rest we have used the relation with other 
measured electron temperatures in an inverse way as de- 
scribed below, 33 from t([Sm]), 20 from t([Oll]) and finally 4 
from t([Nn]). In this latter case, the objects are M51 GEHRs 
of high metallicity from Bresolin et al. (2004). We have as- 
summed as well the t([Om]) to be equal to t([Nem]) and 
t([Anv]) in order to calculate Ne 2+ and Ar 3+ abundances, 
respectively. 

^ IRAF, Image Reduction and Analysis Facility, is distributed 
by the National Optical Astronomical Observatory 



The electron temperature of O , representative of 
the low excitation zone of the ionized gas, has been 
derived from the ratio of the [On] emission lines 
(I(3727A)/I(7319A)+I(7330A)) for a subsample of 311 ob- 
jects. The auroral lines of [On] present higher uncertainties 
because of their lower signal-to-noise ratio and a higher de- 
pendence on reddening correction, due to their larger wave- 
length distance to the closest hydrogen recombination emis- 
sion line. Besides, they present a small contribution due to 
recombination, although in quantities not larger than the 
usual reported errors. Additionally, this ratio has a strong 
dependence on electron density, which makes the determi- 
nation of this temperature and, hence, the determination of 
+ abundances very uncertain. This can increase the un- 
certainty of total oxygen abundances, mostly in high metal- 
licity - low excitation HII regions (Perez-Montero & Dfaz, 
2005). For the rest of the objects, we have used the grid of 
relations between t([On]) and t([Om]) presented in Perez- 
Montero & Dfaz (2003), which takes into account the de- 
pendence of t([On]) on density and its corresponding un- 
certainty. For the same 4 objects of M51 for which neither 
t([Om]) nor t([SIII)]) have been measured we have taken 
t([Nn]) as representative of the electron temperature of the 
low excitation zone using the ratio of [Nil] emission lines 
(I(6548A)+I(6584A))/I(5755A). 

In the case of the electron temperature of Ar 2+ , we 
have assumed that the ionization structure of this ion is 
rather similar to that of S 2+ , whose electron temperature 
has intermediate values between the high and low excita- 
tion zones (Garnett, 1992). We have calculated directly the 
temperature of S 2+ from the ratio of [Sill] emission lines 
(I(9069A)+I(9532A))/I(6312A) for a subsample of 299 ob- 
jects. For the rest of objects we have taken relations of 
this temperature with t([Om]). In the case of low excitation 
objects (T e < 10000 K), based on photoionization models 
(Perez-Montero & Dfaz, 2005): 

t([Sm)] = 1.05 ■ t[{Om)] - 0.08 (1) 

while for the high excitation objects we have taken the em- 
pirical relation found by Hagele et al. (2006) that accounts 
better for the high dispersion found in this relation for HII 
galaxies: 

t([Sm]) = (1.19 ± 0.08) ■ t([Om]) - (0.32 ± 0.10) (2) 
3.2 Ionic abundances 

Each ionic abundance has been calculated using the most 
prominent emission lines for each ion and the appropriate 
electron temperature. In the case of + , we have used the 
following expression using the intensity of the [On] line at 
3727 A, the electron temperature t([0n]), or any other value 
associated to the low excitation zone and the electron den- 
sity, n: 

,o . fO + \ , //(3727)\ r „„„ 1.583 

12 + 109 [w ) = log [im ) + 5 - 992 + — - 

-0.681 -logt + log(l + 0.00023 ■ n) (3) 

In the case of 11 of the 12 very low metallicity HII 
Galaxies identified by Kniazev et al. (2003) and also for a 
subsample of 183 HII Galaxies compiled from Izotov et al. 
(2006), which have no data on the [On] A 3272 A line, the 
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+ abundances have been measured from the emission lines 
of [On] at 7319 and 7330 A using the following expression: 



12 + log 



, / 1(7325) \ „ onr 2.44 

= io 4w) +6 - 895 + T- 

-0.58 • log t - log(l + 0.0047 ■ n) (4) 

For 2+ abundances we can use the [Oni] emission lines 
at 4959 A and 5007 A, along with the electron temperature 
of [Oni]: 



12 + log 



Q 2 + 
H+ 



= log 



[(4959) + 1(5007) \ 

I(H0) ) 

1.251 
+ - 



t 



+ 6.144 + 



0.55- log t (5) 



These two ionic abundances allow to calculate the total 
abundance of oxygen in relation to hydrogen via the expres- 
sion: 

o _ o+ + o 2+ 

H ~ H+ 

Abundances of Ne 2+ have been calculated from the 
[Nem] emission line at 3869 A and the electron tempera- 
ture of [Om] for a subsample of 773 objects: 



(6) 



12 



Ne 2 + 
H+ 



+ _o.504 -log t 



(7) 



Finally, regarding Ar, we can measure Ar 2+ , from the 
[Arm] emission line at 7136 A and the electron temperature 
of [Sm] for a subsample of 572 objects using the following 
expression: 



12 + log 



Ar 2+ \ 



f £(7136) \ 
^ I (Hp) J 
0.808 



+ 6.157 



0.508- log t 



(8) 



It is possible to measure as well the lines of [Anv] at 
4713 and 4740 A. Nevertheless, the first one usually appears 
blended with a line of Hel at 4711 A that is difficult to 
correct, so it is better to use the second one to calculate the 
ionic abundance of Ar 3+ through the electron temperature 
of [Om]: 



12 



M 'Jm\ +4 . 7 „5 + 



HHP) J 
1.246 



0.156- log t 



(9) 



This abundance has been calculated for a subsample of 
253 objects with a simultaneous measurement of the [Arm] 
and [Anv] emission lines. 



4 PHOTOIONIZATION MODELS 

A large grid of photoionization models has been calculated 
in order to check the validity of the ionization correction 
factors for Ne and Ar. 

We have used the photoionization code Cloudy 06.02 
(Ferland et al., 1998), taking as ionizing source the spectral 
energy distributions of O and B stars calculated with the 




0,4 0.6 

o 2+ /(o + +o 2+ ) 

Figure 1. Different relations between the ICF(Ne 2+ ) and the 
ratio of 2+ /0, approximated by 2+ /(0++0 2+ ). The points 
correspond to the models described in the text for stellar atmo- 
spheres of 45 kK and 50 kK. The thick solid line represents the 
best fit to these points. The thin solid lines correspond to differ- 
ent fits of photoionization models as a function of metallicity in 
Izotov et al. (2006). Finally the thick dashed line represents the 
classical approximation. 



code WM-Basic (version 2.11 II, Pauldrach et al., 2001). We 
have assummed spherical geometry, a constant density of 
100 particles per cm 2 and a standard fraction of dust grains 
in the interstellar medium. We have assummed as well that 
the gas has the same metallicity as the ionizing source, cov- 
ering the values O.O5Z , 0.2Z Q , 0.4Z Q , Z Q and 2Z Q , taking 
as the solar metallicity the oxygen abundance measured by 
Allende-Prieto et al. (2001; 12+log(0/H) = 8.69). The rest 
of ionic abundances have been set to their solar proportions. 
A certain amount of depletion has been taken into account 
for the elements C, O, Mg, Si and Fe, considered for the for- 
mation of dust grains in the code. Although it is expected 
that the dust-to-gas ratio scales with metallicity (Shields & 
Kennicutt, 1995), we have checked that different values of 
this ratio do not lead to variations in our results about the 
computed ICFs in quantities larger than the reported errors. 
Regarding other functional parameters we have considered 
different values of the ionization parameter (log U = -3.5, 
-3.0, -2.5 and -2.0) and the stellar effective temperature (T* 
= 35000 K, 40000 K, 45000 K and 50000 K). This gives a 
total of 80 photoionization models to cover the conditions 
of different ionized gas nebulae. Atomic data, including col- 
lision strengths in the models are consistent with those used 
in the calculation of chemical abundances. We have checked 
as well the influence of varying dielectronic recombination 
rate coefficients in the code according with the most recent 
values (Badnell, 2006 and references therein), but we have 
not found any relevant variation in our results concerning 
the calculation of ICFs. 



II Available at http://www.usm.uni- 

muenchen.de / people / adi/Programs /Programs, html 
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2+ + 2+ 

S 7(S +S ) 



Figure 2. Representation of the ICF for Ar as a function of the 2+ /(0 + +0 2+ ) ratio in the two upper panels, and as a function of 
S 2 +/(S + +S 2 +) in the lower ones. Dashed lines represent the expressions given by Izotov et al. (1994), the thin solid lines represent 
those in Izotov et al. (2006). Finally, the points represent the models described here for cluster effective temperatures of 45 kK and 
50 kK. The thick line is the best fit to these points. 



5 DISCUSSION 

5.1 ionization correction factors (ICF) 

ICFs stand for the unseen ionization stages of each element. 

J = ICF(X + 1 ■ (10) 

where X is the element whose ICF is required and X +% 
is the ionic species whose abundance is calculated by means 
of the detected emission lines. 

In the case of neon, only the [Nelll] line at 3869 A is 
detected in the optical spectrum in a large sample of ob- 
jects. Since its ionization structure is quite similar to that 
of oxygen, the following approximation can be used: 



that leads to: 

ICF(Ne 2+ ) = 0/0 2+ « (0 + + 2+ )/0 2+ (12) 



This relation is shown as a thick dashed line in Fig- 
ure 1. Nevertheless, Izotov et al. (2004) point out from a 
set of photoionization models that this approximation is 
less accurate in regions with a lower ionization parameter 
where the charge transfer reaction between 2+ and H° be- 
comes more efficient. In this type of objects it is expected 
to find larger abundances of Ne 2+ in relation to 2+ and, 
hence, the 2+ /0 ratio provides larger ICFs. This trend is 
confirmed with the set of photoionization models used by 
Izotov et al. (2006) that, besides, find a slight dependence 
of the ICF(Ne 2+ ) on metallicity, with larger values of this 
ICF for high metallicity regions. This dependence is due to 
the addition of X-rays sources in their models with lower 
metallicity and not to a real dependence on metallicity or 
ionization parameter. The ICFs from Izotov et al. (2006) are 
shown as thin solid lines in Figure 1. The fit of our models 
for the hottest stars (45000 K and 50000 K), whose effec- 
tive temperatures reproduce better the radiation field in the 
most massive ionizing cluster is shown as a thick solid line 
in Figure 1. It reveals the same overestimate of the classical 
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approximation in low excitation objects. However, we do not 
find any noticeable dependence on the metallicity. This fit 
gives: 



iCF(Ne 2+ ) = 0.753 + 0.142 • x + 0.171/a; 



(13) 



where x = 2+ / (0++0 2+ ), with a rms of 0.074 only. 

Regarding argon, the ICF can vary depending on the 
availibity of Ar 3+ abundances. Izotov et al. (1994) propose 
the following formula to calculate the total abundance of Ar 
using the ionic abundances of these two ions: 



ICF{Ar 2+ + Ar 3+ ) 



0.99 + 0.091 



-1.14 



O 4 



0.077 



(14) 



V o 

'0+\3- - 

O J ' V O ) . 
On the other hand, if only the emission line of [Arm] is 
available, they propose the following expression: 



ICF{Ar 2+ ) = 



15 + 2.39^) -2.64^ 



(15) 



Both formulae are shown as dashed lines in the two 
upper panels of Figure 2, respectively. Nevertheless, these 
fits clearly overestimate the amount of total Ar compared 
with the expressions proposed by Izotov et al. (2006) who, as 
in the case of Ne, find a dependence with metallicity. These 
fits are shown in the upper panels of Figure 2 as thin solid 
lines for low, intermediate and high metallicities. Finally, 
the fits to our own models are closer to these latter, but 
we do not find any relevant dependence of our models on 
metallicity. The expressions for the fits of these models are, 
as a function of x = 2+ / (0 + +0 2+ ): 

ICF (Ar 2+ +Ar 3+ ) = 0.928+0.364-(l-a:)+0.006/(l-a;)(16) 

with a rms of 0.011, shown as a thick solid line in the left 
upper panel of Figure 2, and 

ICF(Ar 2+ ) = 0.596 + 0.967 ■ (1 - x) + 0.077/(1 - x) (17) 

with a rms of 0.067, shown as a thick solid line in the right 
upper panel of Figure 2. It is possible, as well, to express 
these ICFs as a function of sulphur abundances. This allows 
to derive total argon abundances using red or near-IR obser- 
vations only, in the case of ICF of Ar 2+ . Besides, these ICFs 
can be used in the subsample of SDSS objects whose [On] 
line at 3727 A is not detected but which, on the contrary, 
have good detection of the [Sm] at 9069 A. The fits of our 
models, as a function of x = S 2+ / (S + +S 2+ ) are: 

ICF(Ar 2+ +Ar 3+ ) = 0.870+0.695-(l-a;)+0.0086/(l-:r)(18) 

with a rms of 0.019, which is shownd as a thick solid line in 
the left lower panel of Figure 2, and 

ICF(Ar 2+ ) = 0.596 + 0.967 ■ (1 - x) + 0.077/(1 - as) (19) 

with a rms of 0.068, shown as a thick solid line in the right 
lower panel of Figure 2. 



5.2 Behaviour of Ne/O and Ar/O with metallicity 

We have applied the ICF(Ne 2+ ) derived from our set of pho- 
toionization models to a total sample of 578 HII Galaxies, 
117 GEHRs and 12 HII regions in our Galaxy and the Mag- 
ellanic Clouds. In Figure 3, we represent the ratio of Ne/O as 
a function of the total oxygen abundance, and we compare 



0.5 



(),() 



z 

o -i.o 



-1,5 



T T i tTi' 1 t #4-hTT. rLlnTTl 




-2,0 

7,0 7,2 7,4 7,6 7,8 8,0 8,2 8,4 8,6 
12+log(0/H) 

Figure 3. Relation between the abundance ratio Ne/O and the 
metallicity for the sample of objects described in the text. Solid 
circles represent HII galaxies, open squares Giant Extragalactic 
HII regions and open triangles HII regions in the Galaxy and the 
Magellanic Clouds. We show as well the solar chemical abundances 
for O (Allende-Prieto et al. 2001) and for Ne (Grevesse & Sauval, 
1998), linked by a solid line with the solar Ne/O ratio from As- 
plund et al. (2005). The dashed line represents the ratio calculated 
by Drake & Testa (2005) for a set of stellar coronae. 



the obtained values with the solar one, taking the solar oxy- 
gen abundance to be 12+log(0/H) = 8.69 (Allende-Prieto et 
al., 2001) and the solar neon abundance to be 12+log(Ne/H) 
= 8.08 (Grevesse & Sauval, 1998). Although there is a high 
dispersion, probably due to chemical inhomogeneities and 
different observational conditions, the values agree quite well 
with the assumption of a constant value of Ne/O, at least 
for low and intermediate metallicities. At higher metallici- 
ties, there is a slight decrease of Ne/O with O/H, perhaps 
due to the underestimate of the corresponding ICF in this 
regime, as pointed out by Izotov et al. (2006). Neverthe- 
less, since we have not used additional X-ray sources in our 
models, this problem has to be further investigated. At this 
point, it is necessary to stress that all the results from pho- 
toionization models are in contradiction with Vermeij & Van 
der Hulst (2002) who find Ne + abundances even larger that 
those predicted by the classical approximation in a set of 
HII regions in the Galaxy and the Magellanic Clouds from 
direct ISO observations of the emission lines of [Nell] and 
[Nem] in the mid-infrared. Regarding the constant value of 
Ne/O, the average value of the sample, which is -0.72 ± 
0.13, is lower, although within the error, than our assumed 
solar one, but higher than the value reported by Asplund et 
al. (2005) , connected to the value assumed here with a solid 
line in Figure 3 (log (Ne/O) = -0.82 in the photosphere). 
Nevertheless, all the derived values are far from the value 
found by Drake & Testa (2005) using Chandra X-ray spec- 
tra in the coronae of 21 stars, which is log(Ne/0) ft! -0.4. If 
this value is correct, then there is a clear underestimate of 
neon abundances as derived from optical collisional lines. 

Regarding the relation between the Ar/O ratio and oxy- 
gen abundance, our results are shown in Figure 4. In order 
to calculate total argon chemical abundances we have used 
ICFs based on the 2+ /0 ratio for (Ar 2+ + Ar 3+ ) (104 
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Figure 4. Relation between the abundance ratio Ar/O and the 
metallicity for the sample of objects described in the text. We 
show as well the solar chemical abundances for O (Allende-Prieto 
et al., 2001) and for Ar (Grevesse & Sauval, 1998) linked by a 
solid line with the solar Ar/O ratio from Asplund et al. (2005). 
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Figure 5. Relation between the abundance ratio Ar/O and the 
metallicity for the disks of some spiral galaxies. From left to right 
and up to bottom: M101, M51, M33 and NGC 2403. 



objects) and Ar 2+ (344 objects), and ICFs based on the 
S 2+ /(S + +S 2+ ) ratio for the SDSS galaxies from Izotov et 
al.(2006) which do not have observation of the [On] 3727 
A line but for which measurements of the [Sin] 9069 A line 
exist: 44 galaxies in the case of (Ar 2+ + Ar 3+ ) and 101 in 
the case of Ar 2+ . The value of Ar/O presents a larger dis- 
persion than in the case of Ne/O, with an average value sim- 
ilar to the solar value from Allende-Prieto et al. (2001) and 
Grevesse & Sauval (1998), represented by the solar symbol 
in Figure 4 and slightly higher than the value from Asplund 
et al. (2005), connected to it with a solid line. The relation 
appears to be slightly different for HII galaxies, whose be- 
haviour is quite similar to that in the Ne/O diagram, and 
for Giant Extragalactic HII regions, for which there exists a 
slight trend of decreasing Ar/O for increasing metallicites. 



A more careful analysis for individual disk galaxies, as can 
be seen in Figure 5, shows this trend more clearly, contrary 
to it is expected for the production of an alpha element like 
Ar. This is compatible with the existence of larger Ar/O 
at lower effective radius, since there is evidence of negative 
gradients of metallicity in all these galaxies (M101: Ken- 
nicutt & Garnett, 1996; M51: Diaz et al, 1991; M33 and 
NGC2403: Garnett et al., 1997). This result agrees quite 
well with that obtained for sulphur, both for Extragalactic 
HII regions (Diaz et al., 1991; Perez-Montero et al., 2006) 
and halo metal-poor stars (Israelian & Rebolo, 2001). 



5.3 Empirical parameters based on Ne and Ar 
lines 

Although the [Nem] emission line at 3869 A and the [Arm] 
emission line at 7136 A are fainter than the oxygen emission 
lines, which are commonly used in different empirical cali- 
brations of chemical abundances, they can be useful to as- 
certain metallicities when these lines are not available, either 
because they shift out of the optical region due to the object 
redshift or because the instrumental configuration does not 
cover the blue-green region of the spectrum. 

The empirical calibrators are commonly used in objects 
whose low signal-to- noise and/or high metallicities do not 
allow the accurate measurement of any of the auroral lines 
and, therefore, it is not possible to derive the electron tem- 
perature and the ionic chemical abundances from the strong 
collisional lines. 

This is the case for the emission line ratio I([Nem] 
3869A)/I([Oll] 3727A), proposed by Nagao et al. (2006) as 
an empirical calibrator useful for high-redshift galaxies (up 
to z=1.6 in the optical part of the spectrum) and relatively 
independent of reddening due to the proximity of the two 
lines. We show in Figure 6 the relation between this ratio and 
the oxygen abundance for the sample of objects described 
in Section 2. This Figure shows a very high dispersion in 
all the metallicity range but which is especially large in the 
high metallicity regime. This is confirmed by the objects in 
the Galaxy and the Magellanic Clouds compiled by Peim- 
bert et al. (2006) for which it exists a direct determination 
of the oxygen abundance by means of recombination lines. 
The residuals between the direct determination of the oxy- 
gen abundance and those obtained from the relation based 
on the [Nem]/[Oii] ratio are shown in the right panel of the 
same Figure. The standard deviation of these residuals in 
the whole metallicity regime reaches 0.83 dex. 

The reason for this huge dispersion can be found in 
the high dependence of the [Nem]/[0 n] ratio on ionization 
parameter, as it is stressed by the same authors that propose 
it. This dependence is due to the tight relation between the 
flux of [Nem] and [Om] emission lines (see Figure 7), as a 
consequence of the quite similar ionization structure of 2+ 
and Ne 2+ and the constant value of the Ne/O ratio. The 
best linear fit between the fluxes of the lines for this sample 
yields: 

J([Om]4959 + 5007A) = (15.37 ±0.25) • I([7Yein]3869A)(20) 

Although, there is a deviation from this trend for some 
HII galaxies showing abnormally high values of the line of 
[Nem] (e.g. UM382; Terlevich et al., 1991, HS1440+4302, 
HS1347+3811; Popescu & Hopp, 2000), most of the objects 
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Figure 6. Relation between the logarithm of the ratio of the 
[Nem] 3869 A and [On] 3727 A lines and the oxygen abundance 
for the sample described in Section 2, along with a subsample 
of objects in our Galaxy and the Magellanic Clouds with oxygen 
abundance determinations from oxygen recombination lines (solid 
diamonds). In solid line, it is shown the relation proposed by Na- 
gao et al. (2006) and in the right panel the residuals between the 
directly derived abundances and those obtained from this calibra- 
tion. 
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Figure 7. Relation between the intensity of the emission line of 
[Nem] at 3869 A and the sum of the lines of [Om] at 4959 and 
5007 A for the sample of objects described in Section 2. The solid 
line represents the best linear fit to the sample. 



are quite close to this relation. Even for the objects with 
low intensities, which correspond to high metallicity/low ex- 
citation regions and for which the charge transfer reaction 
between 2+ and H° becomes more important, there ex- 
ists a very good agreement with this linear fit. A possible 
cause to the deviation from this relation in some HII galax- 
ies comes from the fact that Ne + has an ionization potential 
5.85 eV larger than + and therefore is more sensitive to 



the stellar effective temperature, which is higher in these HII 
galaxies. One of the consequences of this relation is that the 
[Nelll]/[On] ratio is, in fact, almost equivalent to [Oni]/[Oll] 
which is highly dependent on ionization parameter and effec- 
tive temperature (Perez- Montero & Diaz, 2005) and presents 
a lower dependence on metallicity. 

The relation between [Om] and [Nem] can be used as 
well in all the other diagnostic ratios involving [Om] lines 
in a similar way. This is the case of the O23 parameter, also 
known as R23, defined by Pagel et al. (1979) as the rela- 
tive sum of the strong lines of [On] at 3727 A and [Om] at 
4959 and 5007 A in relation to H/3. The relation between 
this parameter and the oxygen abundance is widely used for 
objects at large redshifts due to the relatively blue wave- 
length and high intensities of the involved lines. There exist 
many different calibrations of this parameter, whose main 
drawbacks are well known: (i) its relation with metallicity 
is double-valued with increasing values of O23 for increas- 
ing values of metallicity in the low metallicity regime and 
decreasing values of O23 for increasing metallicity in the 
high metallicity regime, requiring external methods to dis- 
tinguish the upper from the lower branch and with a high 
dispersion in the middle-range, that makes metallicites quite 
uncertain in the range 8.0 < 12+log(0/H) < 8.4; (ii) the 
dependence of O23 on other functional parameters like ion- 
ization parameter or effective temperature. This is solved 
with the calibration of the parameter as a function of other 
quantities that reduce this dependence, as in the case of the 
[0ll]/[0m] ratio (Kobulnicky et al., 1999 from the models of 
McGaugh, 1991) or the P parameter (Pilyugin, 2000); (iii) 
the lack of objects with directly derived abundances in the 
high metallicity regime makes the upper branch calibration 
to rely heavily on photoionization model results, which dif- 
fer appreciably depending on the model atmospheres used 
(Morisset et al. 2004) and the chosen input conditions. Be- 
sides, these models predict higher metallicities than those 
derived from electron temperatures of [Nil] or [Sm] in inner 
disk regions (e.g. Bresolin, 2006), but not with those derived 
from oxygen recombination lines. The difficulties found to 
derive accurate metallicites using this calibration and, in 
general, all the other empirical relations, make them appro- 
priate to study distributions of metallicity in a statistical 
way in large surveys of emission lines objects, but they do 
not offer reliable results for individual determinations. 

We can then define a proxy for the O23 parameter, us- 
ing the relation between [Om] and [Nem] lines, previously 
discussed: 



02JVe3 = 



I([Qn]3727A) + 15.37 ■ I([iVeiii]3869A) 
I(H(3) 



023(21) 



that is shown in Figure 8 for the sample described in Sec- 
tion 2. We also show as a solid line the equivalent to the 
calibration of the O23 parameter based on the models of Mc- 
Gaugh (1991), for an average value of the ionization param- 
eter (log([Om]/[Oll]) = 1). This calibration gives the lowest 
dispersion in relation with the sample of objects in both the 
lower and the upper branch (Perez- Montero & Diaz, 2005). 
We also show the sample of objects with a direct determina- 
tion of the oxygen abundances based on recombination lines 
in order to illustrate how this parameter has the same prob- 
lems described above for O23. We can redefine the 02iVe3 
parameter relative to the closest and brightest hydrogen re- 
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Figure 10. Residuals between the oxygen abundances derived from the McGaugh (1991) calibration of O23 applied to 02jVe3 
(left) and to 02Ne3' (right) and the oxygen abundances derived from the direct method. Upper panels show the calibration for 
low metallicites and lower panels for high metallicities. 
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Figure 8. Relation between oxygen abundance and the 02iVe3 
parameter for the sample of objects described in the text (sym- 
bols) and the calibration of McGaugh of the O23 parameter for 
an average value of the ionization parameter (log([Om]/[Oll]) = 
1, solid line). 



combination line, H<5. These lines are closer in wavelength, 
which makes its relation to be little reddening and flux cali- 
bration dependent, and useful to larger redshifts (up to z ~ 
1.3). 



O, 



7([On]3727l) + 15.37 ■ J ([/Vein] 3869 A) 



I(HS) 



2 



I {HP) 



(22) 



I(HS) 

This parameter being equivalent to O23, it suffers from 
its same problems, with the additional difficulty of being 
based on weaker lines and therefore, more difficult to mea- 
sure with good signal-to-noise ratio. In the case of the hydro- 
gen recombination line H<5, it is affected more importantly 
by the presence of the absorption line due to underlying stel- 
lar populations and its measurement has to be carried out 
carefully (e.g. Hagele et al., 2006). 
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Figure 9. Relation between oxygen abundance and the C>2jVe3' 
parameter for the sample of objects described in the text (sym- 
bols) and the empirical calibration of O23 of McGaugh for an 
average value of the ionization parameters (log([Oin]/[Oll]) = 1, 
solid line). This calibration has been displaced in the x axis 0.6 
dex to take into the account the mean value of the H/3/H<5 ratio. 



We show this parameter in Figure 9 for the sample of 
objects described in Section 2 and the sample of HII re- 
gions with a determination of the oxygen abundance based 
on recombination lines. We show as well as a solid line the 
empirical calibration from McGaugh but in this case we take 
into account the factor H/3/H5 in the x axis. As in the case 
of O23, it is necessary to distinguish between the upper and 
lower branch to derive the oxygen abundance from it. In the 
case of O23, this is usually done by means of the [Nil] lines, 
but they are not available in the blue part of the spectrum. 
Other possibility is using the [Nem]/[Oll] ratio, although it 
is very uncertain. From Figure 7 we can see that while ob- 
jects with log([Nem]/[Oll]) < -1.0 probably have high metal- 
licities, the contrary is not true. The residuals of the oxygen 
abundances derived with these new parameters as a func- 
tion of the abundances measured using the direct method 
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log(I(3727)/(IH5)) log 2Ne3 , 

Figure 12. Diagnostic diagrams based on neon emission lines for the sample of objects described in Section 2. At left, the relation 
between the ratios [Nem]/H<5 and [Oll]/H<5 and at right the relation between the ratio [Ncin]/[Oll] and the 02jv e 3' parameter. Solid 
lines represent the analytical division between star forming objects and AGNs. Dashed lines represent the limit of the bands of 
uncertainty 0.15 dex both sides of each relation. 




Figure 11. Relation between the Ar3C>3 parameter and the oxy- 
gen abundances derived from the direct method, along with the 
calibration proposed by Stasifiska (2006) for this parameter. In 
the right panel, we represent the residuals between the directly 
derived abundances and those derived from the calibration. 



for both 02jve3 and 02jv e 3' parameters. However, in the low 
metallicity regime, 12+log(0/H) < 8.0, the dispersion of the 
residuals in relation with the direct method are 0.24 dex for 
the 02jve3 parameter and 0.22 dex for the 02jv e 3' param- 
eter, which are much higher than the 0.13 dex calculated 
for the O23 parameter in Perez-Montero & Diaz (2005). In 
this case the residuals between this calibration and those 
based on the neon emission line are 0.14 dex and 0.16 dex, 
respectively. 

Among the other empirical calibrators based on Ne or 
Ar emission lines, we can analyse also the ArsOs parame- 
ter, defined by Stasihska (2006) as the ratio of the fluxes 
of the emission lines of [Arm] at 7136 A and [Om] at 5007 
A. In Fi gure 11, we show the relation between the oxygen 
abundances and the logarithm of this parameter for the sam- 
ple of objects described in Section 2 having a measurement 
of the [Arm] line. The standard deviation of the residuals, 
which are shown in the right panel of the same Figure as a 
function of metallicity, give a value of 0.35 dex. Besides, the 
behaviour of this parameter in the high metallicity regime is 
quite unclear due to the position of some of the objects with 
a direct determination of the oxygen abundance based on 
recombination lines, that show lower values of Ar 3 03 than 
those expected for their metallicity. 



for both the lower and the upper branches are shown in Fig- 
ure 10. Since we can see, as in the case of O23, both the 
upper and lower branch calibrations show large differences 
with the derived metallicities in the turnover region (8.0 < 
12+log(0/H) < 8.4). The standard deviation of the resid- 
uals in the other ranges are quite similar to those obtained 
for the O23 parameter in Perez-Montero & Diaz (2005) for 
the McGaugh calibration of the upper branch (0.20 dex for 
2 jve3 and 0.18 dex for 2JVe 3' in the 12+log(0/H) > 8.4 
range). The standard deviation of the residuals between the 
metallicites derived from the O23 parameter are 0.06 dex 



5.4 Diagnostic methods for the ionizing source in 
emission line objects 

Strong lines are used also in emission line-like objects to 
find out the nature of the main ionization mechanism: pho- 
toionization due to the absorption by the surrounding gas 
of UV photons emitted by massive stars, power law spectral 
energy distribution associated with active galactic nuclei or 
shock excitation. Since some of these diagnostic diagrams 
for the blue part of the spectrum use the emission lines of 
[Om] at 5007 and 4959 A, we could use as well the [Nem] 
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3869 A emission line and the closest and brightest hydrogen 
recombination line to redefine these diagrams that would be 
useful up to larger redshifts and would also be less redden- 
ing dependent. This is the case of the relations proposed by 
Lamareille et al. (2004), as for instance the relation between 
([Om] 5007A/H/3) and ([On] 3727A/H/3). Using the appro- 
priate relations, the analytical expression proposed by these 
authors take the following expression: 



, .' Weill] 3869 A \ 

log I — m — ) 



0.14 



log([On]3727A/#5) - 2.05 



+0.37(23) 



In the left panel of Figure 12 we show this diagram for 
the objects described in Section 2. Starbursts galaxies and 
star forming regions lie below the solid line. This diagnos- 
tic ratio is, in fact, equivalent to the relation proposed by 
Rola et al. (1997) between ([Nelll] 3869A/H/3) and ([On] 
3727A/H/3). The rate of coincidences between this diagnos- 
tic diagram and its predecessor based on [Om] lines reaches 
83 %. This diagram can be used also to separate Seyfert 2 
galaxies from LINERS using the theoretical expression pro- 
posed by Lamareille et al. (2004) adapted for [Nelll] and 
R8: 



log 



[iVein] 



0.04 



(24) 



All the HII galaxies which present abnormally high val- 
ues of the emission line of [Nelll] shown in Figure 7 appear in 
this diagram in the Seyfert 2 zone, what is consistent with a 
possibly higher dependence of the [Nelll] line on stellar effec- 
tive temperature, as it has been previously suggested. This 
is supported by the fact that any of these galaxies are clas- 
sified as Sy2 when other diagnostics diagrams, like [Sll]/Ha 
vs. [OlIl]/H/3 or [Nll]/Ha vs. [Oiii]/H/3 are used. 

The other blue diagnostic ratio proposed by Lamareille 
et al. (2004) is the relation ([Om] 4959+5007 A)/([On] 3727 
A) versus O23. The analytical expressions proposed by these 
authors for this diagram can be expressed in terms of [Nelll] 
and H5 adopting the form: 



/ [jVem]3869A \ L5 

09 V [On]3727A ) ~ log 2Ne3 



'2JVe3' — 2.3 



+ 1.21 



(25) 



The diagram is shown in the right panel of Figure 12 
and, again, starbursts galaxies and star forming regions lie 
below this relation, which is shown as a solid line. Although 
the rate of coincidences between this relation and the one 
based on [On] and [Om] is very high (92 %), a larger number 
of objects is found lying on the AGN region as compared to 
other diagnostics which somewhat questions its application. 



6 SUMMARY AND CONCLUSIONS 

We have performed an analysis of a large sample of emission 
line objects with a direct determination of the oxygen abun- 
dance, via the calculation of the electron temperature, and 
an accurate measurement of the emission lines of [Nelll] at 
3869 A and [Arm] at 7136 A. We have recalculated oxygen, 
neon and argon abundances taking into account the electron 
temperature most representative for each ionic species. The 
total chemical abundances have been calculated with the aid 
of new ionization correction factors for neon and argon based 



on a new grid of photoionization models computed using 
Cloudy v06.02 and WM-Basic stellar model atmospheres. 

The new ICF for Ne yields lower abundances for low 
excitation objects in relation with the classical approxima- 
tion 2+ /0 ~ Nc 2+ /Ne, which does not take into account 
the charge transfer reaction between 2+ and H°. This new 
ICF agrees quite well with the fits proposed by Izotov et 
al. (2006) for the high and intermediate metallicity regime. 
Nevertheless, although we do not find any relevant depen- 
dence of this ICF on metallicity as we have not considered 
different X-ray sources in our models, there is some evi- 
dence of the underestimate of the ICF in the high metallic- 
ity regime. Firstly, the values found by Vermeij & Van der 
Hulst (2002) for a sample of Galactic HII regions from the 
Ne emission lines in the mid-IR point to larger abundances 
of Ne + in these objects. Secondly, the study of the Ne/O ra- 
tio as a function of metallicity shows lower values of Ne/O 
in the high metallicity regime. The average value of this ra- 
tio agrees better with the solar abundances than the values 
recently obtained from X-ray observations in stellar coronae 
(Drake & Testa, 2005). 

Regarding Ar, we have obtained new ICFs for both 
Ar 2+ and Ar 2+ +Ar 3+ quite similar to the values found by 
Izotov et al. (2006) but, again, we do not find, any relevant 
dependence on metallicity. According to these new ICFs, the 
values proposed by Izotov et al. (1994) clearly overestimate 
the total abundances of Ar. We propose as well new ICFs 
based on the ratio S 2+ /(S + +S 2+ ), in order to calculate total 
Ar abundances using the red and far-IR wavelength range 
only. The study of the Ar/O ratio as a function of metallicity 
gives contradictory results for HII galaxies and Giants HII 
regions in spirals disks. For the first ones we find a constant 
value of the Ar/O ratio, in agreement with the expected 
results for the stellar production of this element. However, 
for GEHRs we find evidence for decreasing Ar/O with in- 
creasing metallicity. This result is found as well for the S/O 
ratio both in GEHRs (Diaz et al., 1991; Perez-Montero et 
al., 2006) and halo massive stars of our Galaxy (Israelian 
& Rebolo, 2001). Taking into account that the ionization 
structure of Ar and S are quite similar and that they are 
produced in the same stellar cores it is not surprising that 
their ratios behave in similar ways. 

We have studied some empirical parameters of metal- 
licity based on Ne and Ar emission lines. This is the case 
of the [Nem]/[Oll] ratio proposed by Nagao et al. (2006). 
We have shown that this parameter is indeed much more 
sensitive to ionization parameter and effective temperature 
than to metallicity due to the tight relation existing between 
[Nem] and [Om] emission lines. This is due to the constant 
value of the Ne/O ratio and the identical ionization struc- 
tures of Ne 2+ and 2+ . In fact [Om] can be substituted by 
[Nem] in empirical metallicity calibrations and diagnostic 
diagrams. This offers the possibility of extending the range 
of applicability of these relations to objects with higher red- 
shift minimizing at the same time the effects of reddening 
corrections and flux calibration. We have then defined an 
abundance parameter 02jv e 3' equivalent to the commonly 
used O23 but using the sum of the [On] and [Nem] lines 
relative to H5. Although this parameter has the same prob- 
lems as O23 (double-valued, dependence on U, calibration of 
the upper branch) it constitutes a new tool to derive oxygen 
abundances in large deep optical surveys of galaxies up to 
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high redshifts (« 1.3). Using the same principle, we have de- 
fined new diagnostic methods based only on [On] and [Nelll] , 
similar to those proposed by Rola et al. (1997), that could be 
used in the same surveys to separate star forming galaxies 
from active galactic nuclei. The results obtained from these 
diagrams are very similar to those obtained from the other 
ones in the blue part of the spectrum. 
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